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Abstract 

The infrared structure of SU (2) Yang-Mills theory is studied by means of lattice 
gauge simulations using a new constrained cooling technique. This method reduces 
the action while all Polyakov lines on the lattice remain unchanged. In contrast to 
unconstrained cooling, quark confinement is still intact. A study of the Hessian of 
the Yang-Mills action shows that low action (semi-) classical configurations can be 
achieved, with a characteristic splitting between collective modes and higher mo- 
mentum modes. Besides confinement, the semiclassical configurations also support 
the topological susceptibility and generate spontaneous breakdown of chiral symme- 
try. We show that they possess a cluster structure of locally mainly (anti-) selfdual 
objects. By contrast to an instanton or a meron medium, the topological charge of 
individual clusters is smoothly distributed. 
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1 Introduction 



More than thirty years ago, Callan, Dashen and Gross [H [2] (CDG) attempted to derive 
the non-perturbative properties of QCD from classical solutions of the Euclidean equations 
of motion in Yang-Mills theory. Such solution had been found few years before and called 
pseudoparticles or instantons |3], carrying unit winding number (topological charge). The 
program initiated by CDG was only partly successful: while chiral symmetry breaking 
became immediately understandable, confinement remained unexplained in terms of in- 
stantons. Thus, CDG assigned the property of confinement to singular solutions (merons), 
objects with infinite action |1] and half-integer topological charge. The phase of (con- 
fining) QCD was thought to be realised on top of a Kosterlitz-Thouless transition where 
the dissociation of instantons into merons is probed by a action/entropy trade-off. This 
mechanism, however, was beyond the possibilities of a semiclassical calculation and was 
not worked out further at that time. 

Since lattice simulations have become possible, one of the aims was to learn about the mech- 
anism of confinement from the structure of lattice configurations. The simplest method to 
search for a non-perturbative background is cooling [SIE]- Normal cooling makes visible a 
granular structure of action and topological charge. It appeared natural to interprete the 
lumps as instantons. Using improved actions, even unstable instanton-antiinstanton pairs 
became relatively stabilised and observable under cooling [7]. 

The problem is when cooling should be stopped. Ignoring ad-hoc recipes derived from the 
instanton model itself, self-stopping algorithms are preferable. Restricted cooling [SI [9] 
is such an example. It tolerates a certain violation of the classical lattice equation of 
motions in order to locally stop the relaxation. It turns out that this imposes a certain 
scale (other than the lattice spacing) and leads to a scale-dependent characterisation of the 
vacuum structure. Topological structure with a clear scale dependence is also provided by 
spectral filtering techniques based on the lattice Dirac or Laplacian operator [10 } [TT | [T2 | [13] . 
respectively, which do not modify the underlying gauge field. 

So far, the most diligent approach making use of cooling was a combination of blocking 
and inverse blocking. Here, inverse blocking meant cooling on a fine lattice with a perfect 
action, constrained to preserve the long-range structure kept in form of blocked lattice 
configurations [TH [151 UHl [IZ]- In this form, iterated cycles of blocking and inverse block- 
ing avoid the collapse to the trivial vacuum. It automatically finishes with a nontrivial 
"background" configuration stabilised by the blocked configuration. This numerically very 
expensive procedure was later mapped to simple APE-type smearing [18]. By extrapolat- 
ing back to the "no smearing" limit, the instanton content of SU{2) [20] and SU{3) [21] 
pure gauge theory and full QCD with dynamical fermions [22] has been described. 

There is a caveat, however. The attempted reconstruction of the lattice configurations by 
reduction to their "instanton content" (positions and sizes) does not reproduce confinement 
and the correct hadron spectrum [19]. Kovacs [23] showed that other (low- action) degrees 
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of freedom (torons) must be extracted from the individual lattice configurations as well 
and added to the model configurations in order to mimic the original lattice ensemble in a 
satisfactory way. 

At finite temperature, calorons (finite-T instanton solutions) more general than previously 
known have been found and described by Kraan and van Baal [2H[25] and Lee and Lu [26]). 
These solutions have two advantages over instantons: they are explicitely dependent on the 
embedding holonomy, and the moduli space contains the necessary degrees of freedom to 
permit dissociation into Nc constituents, dyons carrying fractional topological charge. One 
of the dyons carries the zero mode that would otherwise assigned to the whole caloron. A 
gas of moderately dissociated calorons is able to predict confinement provided the holonomy 
parameter is maximally nontrivial (i.e. the external Polyakov loop vanishes) [271 128] . 
Although this semiclassical model has not been completely worked out to the extent that 
the effective potential of the Polyakov loop and the confining force could be obtained 
as functions of the temperature, the expected caloron/dyon structures could be actually 
detected on realistic lattice configurations below the deconfinement transition [29 1 l30 l l3Tj. 
In later studies cooling has not been used. Instead, modes of chiral (overlap) fermions were 
used as probes. In the deconfined phase the symmetry between dyons is broken, and only 
the light ones are abundant [32j. 

In a model-independent way, the question of confinement facilitated by semiclassical con- 
figurations has been asked again and affirmatively answered by Gonzalez Arroyo and col- 
laborators [331 1311 I3S]- They have used cooling, however employing twisted boundary 
conditions [361 EH EH]. Under these circumstances a trivial vacuum is avoided and frac- 
tional topological charges can be observed. 

At present, a semiclassical picture for confinement seems not very popular. As far as 
a mechanism of confinement is concerned, the overwhelming opinion is that exclusively 
defects like vortices and monopoles [391 HO] play this role Q. Practically, they are exposed 
as defects of gauge fixing. It is not excluded, however, that certain non-Abelian extended 
objects (e.g. "thick vortices" [12]) exist, which become merely localisedhj the process of 
gauge fixing and projection to the Abelian degrees of freedom. We should mention that 
also calorons lead to a characteristic monopole and vortex structure [2S] . 

A common feature of the last examples, beginning with inverse blocking, is the interplay of 
topological configurations and certain constraints or boundary conditions imposed to the 
cooling process. We think it is worth to pursue this guiding idea further. In a preceding 
paper, one of us has explored the feasibility of constrained cooling that preserves the 
property of confinement [13] by keeping fixed Wilson loop ratios in the given configuration 
at a given "macroscopic" (infrared) length scale. The cooling constrained by large-scale 
(infrared) observables is formulated by adding Lagrange multiplier forces to the lattice 
equations of motion. This constrained cooling has also led to a gas formed by clusters of 
action. Not surprisingly, the cooled configurations reproduce the conventional confinement 

^We refer to Ref. [41] for a way to reconcile the two points of view. 
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criteria despite the low action and the apparent absence of strong fluctuations. 

In this paper we explore the possibilities of another type of constrained cooling. It uses 
the information stored in the (locally varying) Polyakov loop. Because we study first the 
situation at T = 0, we preserve the Polyakov loops in all (temporal and spatial) directions. 

The paper is organised as follows. In the following Sect. 12.11 we consider in general 
terms the connection between classical configurations, constrained cooling and streamline 
configurations. In Sect. 12.21 we explain how cooling constrained by Polyakov loops can 
be realised. Here we report about the easiest available observations, the behaviour of 
action and Creutz ratios under cooling, and contrast this with standard unconstrained 
cooling: confinement is preserved. The next Sect. I2.3l is devoted to the spectrum of gluonic 
excitation modes on the background of streamline configurations and the difference to the 
spectrum on standard-cooled configurations. In Sect. [3] we describe the lumpy structure 
in terms of action and topological charge. In addition, the degree of ( ant i-) self duality and 
the locahsation and dimensionality of action and topological density is analysed. Sect. H] 
describes how the spectrum of the staggered Dirac operator depends on constrained cooling 
and gives a comparison with standard cooling. In the case of constrained cooling we find 
that a set on near- zero modes survives, indicating that chiral symmetry breaking survives. 
In Sect. [5] we draw conclusions for further work. In an Appendix more details about solving 
the constraints during cooling are worked out. 

2 Confinement from Yang-Mills streamlines 

2.1 Classical configurations, streamlines and cooling 

As an illustration, let us consider a scalar field (f){x) with the Euclidean action S[(j)]. Let 
the action possess several degenerate absolute minima which might arise from a discrete 
symmetry of the action. These minima are separated by potential wells over the multidi- 
mensional configuration space. The well known double-well problem [H] might serve as an 
example: there, the two minima of the action are related by the reflection — )■ —0 of the 
corresponding scalar fields. If this symmetry is not spontaneously (or anomalously) bro- 
ken, a semi-classical calculation of the partition function starts by taking into account the 
contributions of all absolute minima weighted by the contribution from Gaussian fluctua- 
tions around each ground state configuration. The implicit assumption for a semi-classical 
treatment is that only configurations close to these absolute vacuum states are relevant 
and tunnelling is rare. It is evident that this treatment breaks down if the potential wells 
between the minima are small and tunnelling becomes frequent. Formally, this can be seen 
that space-time dependent field configurations (e.g. kinks) which interpolate in configu- 
ration space between two possible vacuum states possess an action close to the vacuum 
action implying that, depending on entropy, those kink configurations are relevant as well 
in the full quantum theory. 
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Low action solutions are generically obtained by the so-called cooling procedures which 
minimises the action. If we consider a gas of as many kink as anti-kinks, excessive cooling 
results in an annihilation of kink anti-kink pairs leading to the trivial vacuum. If interesting 
non-perturbative phenomenons can be attributed to the kink configurations, cooling must 
be stopped by either an ad hoc criterion or by additional physical input. 

Reducing the action by means of the steepest descent method, in other words standard 
cooling, can be considered as a flow in configuration space where the field picks up an 
additional dependence on the "cooling time" r, i.e., — t- (/)(r, x). By constrained cooling, 
we address a method which reduces the action while preserving a constraint, such as 

= constant, ^ n^{x) = 1. (1) 

X X 

Both cooling methods can be formally defined by the flow equation: 
It is easy to show that 



d 

n{x) (f){r, x) = 



implying that the constraint has been implemented exactly. Note that the flow equation 
for standard cooling is obtained by setting n(x) = 0. 

Consider (f)^( low action solution for the case the constraint has been switched off, 

i.e., n{x) = 0. If (f)'^{x) is, according to some norm, close to a minimum configuration 
(pmin{x), the difference between (j)'^{x) and 0(r, x), 

A0(r, x) := (f){r,x) — (j)'^{x) , 

is small for some large value r. Note that, with the constraint in place, the constraint 
has to be sufficiently weak such that A0 can be small. Whether this is indeed the case 
for a particular theory has to be justified by a numerical experiment. With the above 
assumption, we might then approximate the action by a second order polynomial: 

S[<P] = S[<p^] + 5^ A0(x)/[0^(x) (3) 

X 

+ IYI ^'^(^) Mm{x,y) A<j>{y) + 0{A<p') . 

xy 

The drift force / and Hessian M are thereby given by 

d(f){x) d(f){x) 0(j)[x) (t>=<j>'' 
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Inserting (|3]) in the flow equation ([2]), we find: 

c?0(r, x) 



dr 



fx - M^y A(l)y + (^Uylfy + My^Acj)^]^ 71^ 



(5) 



where we have introduced the convention to sum over indices which appear twice. Intro- 
ducing the projector 



the last equation can be written in compact notation: 

90(r, x) 



dr 



-Fx - MxyA(j)y , 

p/ , M = PM . 



(6) 
(7) 



The flow equation can be easily solved 

0(r) = 0^= + exp|-A^ r [> 0(0) - l-exp<{-A^r 

where 



G 



MG = F => G = {PMPy^F 

given that PMP is invertible in the projected space. 
We firstly study standard cooling where we have: 

n{x) = 0, M = M = M^, F = f . 

Since the Hessian M is symmetric, we can introduce its eigenmodes as a complete basis, 

Mxy = A„ iljf^ , Ai < A2 < . . . (9) 



I.e. 



and expand 



A0. = J]c„(r)^(") 



(10) 



The solution (JHl) of the flow equation now translates into 



fn 



Cnir) = Cn{0) exp{-A„r} - — 1 - exp{-A„r} 



:ii) 



In the limit r oo, we observe c„ — )■ —fn/^n and approaches the true minimum of the 
action ([3]), i.e., — 0min(a;): 



lim 0(r) 

T—^OO 



(12) 



The crucial idea is to stop cooling by only admitting r < Tstop- As explained by Huang 
in |13], the success of a semi-classical expansion critically depends on the existence of a 



gap in the eigenmode spectrum of the Hessian. Assume that such a gap occurs at the ith 
eigenvalue, i.e., Aj ^ Aj+i. In this case, we will find a range R of t values with 



|r G [Ti,Th] XiTh < 1, \i+lTi > l| . 



(13) 



R = 

For T & R, the higher eigenmodes have died out, 

Cn{T) ~ for n > i + 1 , 

while the lower eigenmodes are "drifting" along the gradients of the action: 

c„(r) = - fnT for n <i . 

Hence, by choosing rgtop G R f|T3|) . it is possible to separate Gaussian fluctuations from 
the low-action collective modes. Those configurations, which are not quite minima of the 
actions, are called streamlines throughout this paper. 

Although it might be feasible in practice to choose Tstop G R hence generating a streamline 
configuration, the properties of those streamline configurations might still depend on the 
ad hoc chosen Tstop- It would be highly desirable to replace the ad hoc stopping of cooling 
by a physical criterion thereby linking the properties of the streamlines to a feature of 
the theory. Constrained cooling is precisely delivering that. To see this, let us return to 
the flow equation ([8]) now with the constraint in place, i.e., n{x) not all zero. We do not 
artificially stop cooling and consider the limit r — >■ oo (although we admit that this is likely 
to be impracticable in an actual computer simulation): 



lim0(r) = 0^ - {PMPy'Pf ^ (j)rain. (M) 
r— >oo 

Note that if the constraint is sufficiently weak in the sense 

{PMPy^Pf ^ M-^ f , 

the limiting configuration would correspond to a streamline configuration of the un-con- 
strained theory. 

Below, we will detail that the above assumptions are met for the case of an SU(2) gauge 
theory where the constraint is in place to conserve the quark confinement. 



2.2 Constrained cooling in Yang-Mills theory 

In Yang-Mills theories, minimal action solutions are the celebrated instanton configurations 
which are characterised according to their integer- valued topological charge Q. Most rele- 
vant at quantum level is the sector with vanishing topological charge: to this sector, a gas 
or liquid of as many instantons (Q = 1) as anti-instantons {Q = —1) would contribute be- 
sides the empty perturbative vacuum. Highly relevant properties of the Yang-Mills vacuum 
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such as the spontaneous breaking of chiral symmetry can be attributed to the instanton 
medium. It is well known that excessive cooling leads to an annihilation of instanton 
anti-instanton pairs or even to charge changing (anti-) instanton elimination. Since we are 
eventually left with the empty vacuum, cooling has to be stopped to retain viable non- 
perturbative vacuum structure. Here, we will show that the constrained cooling detailed 
below results in Yang-Mills streamline configurations which not only give rise to sponta- 
neous chiral symmetry, but also to quark confinement. The stucture, however, differs from 
a gas or liquid of instantons. 

For this purpose, we extract streamline configurations, which are non-perturbatively rele- 
vant, from Monte Carlo configurations occurring within the lattice formulation. Starting 
from a snapshot taken from the Markoff chain simulation of the path integral, we pro- 
pose a constrained minimisation of the action. We shall describe the procedure for SU (2) 
Yang- Mills theory and characterise the emerging field configurations. In the present work 
we concentrate on the zero temperature situation (vacuum state) and use a symmetric 
Euclidean lattice with periodic boundary conditions and finite lattice spacing a. As usual 
for SU (2) gauge theory, the gluon degrees of freedom are represented by unitary matrices 

U^ix) = ulix)+iu^ix)r, ulixf + u^ix)^ = 1, (15) 

associated with the links of the lattice. We are going to work with the Wilson action 

'S'wiison = 5Z ~ ^trP^^(a;) , (16) 

PA^) = U^{x)U,{x + iJi)Ul{x + v)Ul{x) . (17) 
The Yang-Mills partition function is given as the discrete functional integral 

Z ^ jvu exp{-^5wiison}, (18) 

where /3 = 4/(7^ is related to the bare gauge coupling, and represents the relative impor- 
tance of gauge field configurations. This fixes also the length scale by the corresponding 
lattice spacing a = f{(3). The integration measure comprises the Haar measure for each 
hnk. Observables A{U) are calculated by means of a Markoff chain Monte-Carlo approach 
which allows to obtain expectation values by importance sampling of relevant lattice con- 
figurations {U^} and averaging over A{U). 

Starting from such an importance-sampled lattice configuration, a configuration with re- 
duced action is obtained by a mapping summarising the cooling procedure: 

Thereby, we are concentrating on the streamline regime when the Gaussian fiuctuations 
have been removed already. In practice, we will monitor the evolution of observables. 
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averaged over an ensemble of Monte Carlo configurations, as a function of the number 
-^cooi of cooling sweeps. It is well known (and will be demonstrated again below), that 
if cooling is performed by steepest descent with respect to the gauge action, confinement 
fades away with Ai'cooi oo- 

Using the constrained cooling procedure proposed in this paper, our aim is to generate 
streamline configurations from the Monte Carlo configurations, such that an ensemble of 
them, hopefully, still supports the properties of quark confinement and chiral symmetry 
breaking. To this aim, we take a radical attitude as announced above and design a cooling 
procedure which largely reduces the action while it does not allow to change every single 
Polyakov line in all Euclidean directions on the lattice. If P4 denotes one Polyakov line in 
time-like direction, X4 = t, 

Piix) = tr JJC/4(t,f), x^{t,x) 

t 

the static, colour averaged quark-antiquark potential V{r) can be obtained via the corre- 
lator of two Polyakov lines 

exp{-l/(r)7Va} = {P^{x) P^{y)) , \x-y\^r, (19) 

where the average (...) is over an ensemble of configurations {[/}. If (. . .)c denotes the 
average over cooled configurations {V^}, the constraint trivially implies that the cooled 
configurations reproduce the static colour averaged potential 

{P^{x)P,{y))c = {P^{^)PM) = exp{-y(r)7Va} , (20) 

that could also be measured on the ensemble of non-cooled configurations. We will point 
out that not only the static quark-antiquark potential is preserved, but also the poten- 
tial between charges of higher representations. At zero temperature, the four Euclidean 
directions should be equivalent. Therefore we formulate the constraint for cooling such 
that also spatial Polyakov loops are not allowed to change. In the result, one can expect 
that this constrained cooling reproduces also the area law for fundamental Wilson loops 
that gives the colour singlet quark-antiquark potential. This will be checked below in this 
subsection. 

The only way that our original hope could fail is when the constrained cooling fails to reduce 
the action by a significant amount. If the cooling sweeps would be seriously inhibited by the 
constraints, the action would be reduced only by an insignificant amount. This would imply 
that the "streamline configurations" we are looking for would not be smooth semiclassical 
fields. Whether this is the case or not can be decided by a numerical experiment. Below, 
we will show that the specific constraints do not prevent the cooling procedure from largely 
reducing the action. 

The constrained cooling is easy to formulate. Defining the staples as usual by 

= E + I") ^^(^ + ^) ^^^(^) + ^^(^ + - ^) ^^(^ - ^) ^'^(^ - ^) ] ' (21) 



9 




cooling step cooling step 

Figure 1: The action density s (left) and some Creutz ratios Xnn (right) as a function of the 
number of (constrained) cooling sweeps for an ensemble of 784 Monte Carlo configurations 
created at /3 = 2.4 on a lattice 16^. For comparison, the case of standard cooling is also 
shown for s and Xss (1008 configurations). 



the part s^{x) to which a hnk U^{x) contributes is given by 

s^{x) = -ti\u,Xx)B^{x)'^ . (22) 
In order to implement the constraints, we define incomplete Polyakov lines in any direction 

N-l 

'Pf.i^) = l[U^{x + ne^) , (23) 

n=l 

such that the trace of the Polyakov line through some point x can be written as 

P^(f) = tr|[/^(a;)P^(x)}, Vx^ . (24) 

In order to preserve the (hypercubic) isotropy we impose the condition that Polyakov lines 
remain unchanged to all directions fi. The update is done sequentially. If U^^ is the 
link before the update into uIj^^^\ and vji^\x) is the incomplete Polyakov loop with all 
previous updates performed, then the following is required to hold: 

tT[u,{x)V,{x)] = tr{f/J")(a;)p(")(x)} = tr{f;J"+i)(x) P^")(a;)}, Vx, , (25) 

where U^\x) = U^{x). Thus, the original Polyakov loops are preserved. The latter con- 
dition is built in as a constraint using Lagrange multipliers A^(x). When a particular link 
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U^{x) becomes subject to relaxation uj^'' — > W^'^^' , we are led to minimise the principal 



(n+l) 



action 



(26) 

with respect to the link U^{x), the result to be stored as part of the iterated configuration 
(n+i) rpj_^jg .g pai-^ic^iarly easily done for the gauge group SU{2). Thereby B^^\x) is 



the staple formed out of links with all previous updates performed. We find the global 
minimum of the principal action for any choice of the Lagrange multiplier to be attained 
for 



f/("+i)( 



1 

Jf 



4")t(3,) + ^^(^)p(n)t(3.) ^ Ar:=deti/2 B'^^\x) + \^{x)V^;\x) 



(27) 



The remaining task is to fix the Lagrange multiplier such that the constraint in fl25|) is 
satisfied. Details of the calculation are left to the appendix |Al It turns out that there is a 
unique solution for A^. The final result is 



A, 



-^tr(EWtp<")) + ltr(f/J")p(")) AT 



det(4"0 - \i^\B^V, 



(28) 
(29) 



or m summary 



f/i"+i)(x) 



1 



+ 



(30) 



Here no summation over /i is implied. This result has an intuitive interpretation: while 
the term in the square brackets minimises the action locally, the remaining terms 

implement the constraints by projection. 

The constrained cooling method works as follows: 

• generate a thermalised configuration of links {U^j} = {U^^^}] 

• randomly choose a link, say U^{x) on the lattice and "update" the link accordingly 



to local cooling: ujr\x) — )■ Uf 



X] 



if this is repeated 4iV^ times (equal to the 



number of links on the lattice !) this is considered as one cooling sweep; 

• for a sufficiently large (but for practical reasons finite) number A^^cooi, we obtain the 
desired semi-classical configuration. 

In practice, the quantities calculated with these configurations will be carefully monitored 
for a possible residual dependence on A^'cooi- 
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We finally mention a technical point: excessive accumulation of rounding errors can dete- 
riorate the Polyakov line values in the chain 



tr 



tr{f7i")(x)p(^')(x)} ^ tr{[/W(x)pW(x) 



of equalities. For Ai'cooi bigger than (9(10"^), we would advise to modify the coohng step 
iOD by 



M L 



ltr(f/fp(°))pH^ (31) 



det(i?;^"^) - |tr2(i?r'^r) 



\ 1 - itr^(f/rP. 



It is easy to check that the constraints are exactly satisfied: 



Below, we will compare the outcome of our new cooling method with results obtained with 
standard (un-constrained) cooling. The standard local cooling step is given by applying 

f/^"+i)(a;) = Ar:=deti/24")(x) , (standard) (32) 

to a randomly selected link. Again we call 4A^'^ local cooling steps one cooling sweep. 
The configuration reached after A^^cooi sweeps is called "cooled configuration" . We label an 
ensemble of such configurations by the number of sweeps, Ncooi 

Ul{x) = f/f--)(x) , (33) 

even though the actual number of updates some particular link has experienced may differ 
because of the random selection. 

We have tested the constrained cooling approach using an equilibrium ensemble of ther- 
malised configurations generated on a A^^ = 16'^ lattice with Wilson's one-plaquette action 
at /3 = 2.4. For this inverse coupling the lattice spacing a in units of the string tension a 
is roughly given by aa^ = 0.0738(5) [15]. The left panel of Figure [1] shows the ensemble 
average of the space-time averaged action density 

^=]^EE(l-^t^^-(^)) ' (34) 

where P^u{x) is the plaquette ( |T7|) . as a function of the cooling sweep for our constrained 
cooling method and for the standard one. It turns out that, despite of the constraint, 
average plaquette values very close to 1 are reached. We also mention that if the Polyakov 
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lines of a particular configuration are not constant over the orthogonal 3-volume, con- 
strained cooling can never settle with the empty vacuum since this state is incompatible 
with varying Polyakov lines. Indeed, if we fit the action density averaged over the ensem- 
ble of constrained cooled configurations as a function of the number n = A^^cooi of cooling 
sweeps, we find empirically 

s ^ 0.0012(1) + 0.074(1) n-"-'^"^^) , = 0.02 . (35) 

The fit is also shown in figure [T] and seems to be a fair representation of the data for 
n > 100. 

The circumstance that the Polyakov line correlator from the constrained-cooled configura- 
tions equals the correlator from the original Monte Carlo ensemble, implies that the cooled 
configurations share the same quark-antiquark potential with the uncooled SU{2) theory. 
It is also instructive to investigate the amount of disorder remaining in the Wilson loops. 
For this purpose, we consider the Creutz ratios 

(W{t,s))(wit-l,s-l)) 
Xts = - In ) ^^Y^ f , (36) 



where W{t, s) are rectangular Wilson loops with extension t in time and extension s in 
spatial direction. We consider here only diagonal Creutz ratios Xnn- For n such large that 
the excited states have died out, the Wilson loop expectation factorises into an exponential 
term with circumference and area term and a t-independent overlap matrix element. In 
this case, the Creutz ratios exhibit only the area term and approach the asymptotic force 
from below function of n 

lim Xnn = cra^ , (37) 

n—>-oo 

where aa^ is the string tension in units of the lattice spacing. The right panel of Figure [1] 
shows some Creutz ratios Xnn as functions of the cooling sweep. We observe that Xnn 
approaches the string tension for an increasing number n, as it must be. In fact, we observe 
for constrained cooling that the string tension is approached from below with increasing 
n. Without any cooling the Creutz ratios approach the string tension from above. The 
reason for this apparent discrepance is easily understood by noting that the Creutz ratio 
can be interpreted as the derivative of the static quark potential. 



Xnn ~ — 

ar 



1 2 

— + aa . 



Since cooling efficiently wipes out the Coulombic part of the quark potential, there is no 
reason anymore that Creutz ratios tend to overestimate the string tension. For Xnn at 
n = 8, we still observe a slight deviation from the un-cooled value. We attribute this 
observation to a decrease of the gap between between ground state and excited state 
energies during cooling. In this case, even larger values of n (and, in fact, larger lattices) 
would be neccessary to recover the asymptotic value. 
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Figure 2: Left: the spectral density of the Hessian as a function of its eigenvalues for the 
perturbative vacuum and the average spectrum from un-cooled configurations. Right: the 
spectral density after 500 and 1000 sweeps of constrained cooling. Results are shown for 
320 configurations of a 4^ lattice at /3 = 2.3 and one configuration corresponding to /3 = oo 
(perturbative vacuum). 



Standard cooling yields a quite different picture. There Xss is rapidly decreasing and 
the asymptotic value is compatible with zero. This confirms the common wisdom that 
configurations which arise from standard cooling do not sustain quark confinement. 

2.3 How good is the semi- classical expansion? 

Since we have discovered above low action constrained solutions which support important 
features of Yang-Mills theory such as quark confinement, the question arises whether a 
semi-classical expansion around these configurations could be a reliable tool to describe the 
theory even quantitatively. As detailed by Huang in [33] and explained in subsection 12.11 
a first indication would be finding a gap in the eigenvalue spectrum of the Hessian. In 
this subsection, we will show that a clear separation between Gaussian fluctuations (in 
the bulk) and low lying collective modes is clearly observed for our streamline background 
configurations. 

In the context of gauge theories, the parameterisation of the fluctuations must be carefully 
chosen in order to obtain eigenvalues of the Hessian which are invariant with respect to 
gauge transformations of the background field. A suitable parameterisation was given by 
Luescher and Weisz [4B] and will be adopted here: if {f/^} denotes the streamline config- 
uration, i.e., the background field, in the neighbourhood a particular link is decomposed 
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Figure 3: Left: the spectral density of the Hessian as a function of its eigenvalues after 
10000 sweeps of constrained cooling. Right: the spectral density after the same number of 
standard cooling sweeps. The densities are the result of averaging over 320 independent 
4^ configurations. 

into 

U,ix) = exp{tq,ix)}U;ix) , q,{x) = g;(x)r'^ , (38) 

where the q'^{x) represent the fluctuations. These fields change under a gauge transforma- 
tion Q{x) as follows: 

ul^ix) = n{x) u^{x) n\x + fi) , 

Uf{x) = Q{x)U;:{x)Q^{x + f,), (39) 
q^{x) = n{x) q,{x) n\x) . (40) 

Gauge transformations simply rotate the fluctuation field q^{x): 

g;^(x) = R^\x)q'^{x), 
R'\x)t^ = Q{x) t"" Q\x) . (41) 

Expanding the Yang-Mills action up to second order in the fluctuations yields: 

s[u] = s[u^] + $^g;(x)/;r] + ^ E €(^)Ki^^y)9tiy) + o{q'),m 

x,a^ xy,ab,fii' 

where /° is the drift force and M^^ is the Hessian. Because the action is gauge invariant and 
in view of the transformation properties (HOj) . one easily finds that the Hessian transforms 
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homogeneously under a gauge transformation of the background field: 



M;l{x,y)[U^''] = R^%x) M^^(x,i/)[f/1 (R'Y'iy) ■ (43) 

This implies that the eigenvalues of the Hessian are gauge invariant, and the corresponding 
eigenmodes are a mere rotation of the eigenvectors of M^Ulf^]. 

For the Wilson action ( |T6|) and with the definition 

the Hessian is given by: 

MfA^^y) = \ti{u,{x)B,{x)]5^n,,y5,,. 

- \ {^""U^ix) + /i) Ul{x + v) f/t(x)} (5,,,+^ [1 - 5^,,] 

- \ {^"U^[x) U,[x + ^) + v) Ul{x)^ 5,,,+, 5^,, 

- \ {<^"f^M(^) + Vl{x + z/) b''Ul{x)^ 6y,, [1 - 5^,,] 

- I {s^'U^ix) S'Ulix + fi-u) Ulix - u) U,{x - J/)} [1 - 5^,,] 

- I {s'^U^ix) Ulix + fi-u) 6'ulix - v) U,{x - z/)} 5,^, 

- i {5"f/^(x) f/t(x + /i - z/) f/i(x - v) 5'U,{x - z/)} [1 - 5^,,] . 

(44) 

The spectral density p(A) is defined by the number p(A) d\ of eigenvalues A of the Hessian 
in the interval between A and A + c/A. We have investigated this spectral density over 
the whole range of eigenvalues for a small lattice of size since the gross features of the 
spectrum can be already illustrated using rather small volumes. 

Let us firstly study the empty, i.e. perturbative, vacuum where all links are set to unity, 
Ufj_{x) = 1, Va;, fi. The result is shown by the open bars in the left panel of Figure [2l The 
spectrum resembles that of the Laplacian operator besides of a rather large number of zero 
modes. This shows that the empty vacuum possesses a huge number of flat directions. 
This has been previously pointed out in [Ul HHl SH]- In order to imagine one of these 
flat directions, let us assume that all time-like links pointing forward from a given time 
slice is set to an arbitrary SU{2) element, let us say U. In this case, all plaquettes on 
the lattice are still equal to 1. In other words, this configuration is degenerate with what 
we have called perturbative vacuum. On the other hand, the Polyakov line of the new 
configuration is given by tr [/, which is different from tr 1 of the perturbative vacuum. 
Since the Polyakov line is gauge invariant, the new configuration is not gauge equivalent to 
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our starting configuration. Variation of tlie action with respect to tlie degrees of freedom 
of U hence will produce genuine zero modes of the Hessian. 

In the left panel of Figure [21 the full bars show the spectral density averaged over 320 
independent configurations from a thermalised 4^ lattice obtained at /3 = 2.3. The dis- 
crete spectrum from the perturbative vacuum is now washed out. In addition, negative 
eigenvalues naturally appear since a thermalised configurations does not represent a local 
minimum of the action. 

The right panel of Figure [2] show the spectral densities corresponding to different numbers 
of constrained cooling sweeps. Note that there is very little change if the number of cooling 
sweeps is increased from 500 to 1000. As expected, the modes which are present for a 
empty vacuum background get emphasised during cooling. More importantly, there is 
a certain amount of near zero, i.e. collective modes which seem to remain stable under 
constrained cooling. Among those there are also negative eigenvalues of the Hessian. Those 
modes would normally exponentially grow, destroying the streamline configuration. It is 
this growth which is inhibited by the constraints imposed upon the cooling process. 

Note that there is also a gap between the collective low-lying modes and the higher modes. 
The latter we therefore identify with the fluctuations, while the collective modes represent 
collective deformations of the streamline configuration. It is interesting to compare the 
spectrum obtained in the result of constrained cooling with that calculated for configu- 
rations emerging under standard cooling. Our spectral results after 10000 cooling sweeps 
in either case are compared in the left and right panels of Figure [31 The inlays show 
the eigenvalue spectrum near the peak at zero in finer resolution (implying a change in 
the scale due to the fixed normalisation of the spectral distribution). In sharp contrast 
to constrained cooling, the collective modes disappear after a sufficiently large amount of 
standard cooling. 

The gap between the collective near-zero modes and the Gaussian fluctuations at intermedi- 
ate eigenvalues in the bulk spectrum is a generic feature of constrained cooling. According 
to the inspiring paper by Huang [H] , both the presence of near-zero modes and the spectral 
gap between them and the bulk spectrum are the hallmark of semiclassical configurations. 
Admittedly, the SU{2) gauge theory with (3 = 2.3 using an 4^ lattice might be close to 
the finite-volume pseudo-deconfinement transition. We therefore checked the spectrum of 
the Hessian for the same P value using an 8^ lattice. This setting comfortably belongs 
to the confinement phase of the theory. In order to obtain all eigenvalues of the Hessian, 
the matrix has been exactly diagonahsed using the 256 nodes of the dedicated cluster at 
the HPCC, Plymouth. Already the eigenvalue spectrum from a single lattice configura- 
tion shows the familiar pattern (see figure [1]): a clear separation between collective and 
Gaussian modes. 
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Figure 4: The spectral density of the Hessian as a function of its eigenvalues for a single 
configuration prepared at /3 = 2.3 on a 8^ lattice, after 1000 sweeps of constrained cooling. 
The insert shows the vicinity of zero. 



3 Space-time texture 

3.1 Action and topological charge 

Since standard cooling and its refinements (improved action, restricted cooling) have been 
developed with the original intention to analyse deeper the vacuum texture, an inhomo- 
geneous structure visible in action and topological charge, it is clear that we also have to 
further analyse the space-time structure which characterises the configuration after con- 
strained cooling. A particularly important questions is whether the emerging streamline 
configurations are related in any way to instantons, the apparent structures believed to be 
seen in some stage of standard cooling. 

To define the notation, we briefly review the basic properties of an 1-instanton configura- 
tion. In the continuum formulation, the gluon field and the field strength tensor are given 
by (in the regular gauge) 

A;{x) = ^.x'^^^, (45) 



with tJq- = — r^^Q = Sai and rj"-/^ = eaik ■ The action of the instanton is independent of the 



18 



180 I < 1 r 




Figure 5: Left: the field strength defined as the clover average over plaquettes touching in 
the point x: C^u{x) = (Pi + P2 + P3 + Pa)/4:. Right: the topological susceptibility Xtop as 
a function of the cooling sweep for constrained and standard cooling. The configurations 
have been prepared on a 16^ lattice at /9 = 2.4. 



instanton radius p and is given by 

^l-instant = E / = ^ " (^7) 

The (anti-)instanton is an ( ant i-) self dual solution of the Yang-Mills equation of motion: 
F;,(x) = ±F;,(x), = ^e,,^^ F^^ix) . (48) 

We also introduce for convenience the colour electric and colour magnetic fields {i,k,l G 
{1,2,3}): 

F,^^{x) = E^{x) , Ft,{x) = e,uBt{x) . (49) 

Smooth configurations composed out of (anti-)instantons are characterised by their topo- 
logical charge, expressible in terms of the field strength by 

^ = ^ E E/ d'xe,.a^F^^,ix)F:,ix) = J d'x q{x) , (50) 
with the topological density 

^(^) = i E ■ (51) 

i=l,3;a 

19 



lump number 


-^lump 


Qlump 


A, 

^lump 


lump number 


-^lump 


Qlump 


A, 

^lump 


1 


4251 


0.6318327E+00 


0.3244866E+01 


1 


1600 


-0.1349146E+01 


0.1967493E+01 


2 


5123 


0.1609389E+00 


0.1796208E+01 


2 


1536 


0.7005824E+00 


0.1241354E+01 


3 


509 


0.5354339E+00 


0.1027845E+01 


3 


130 


-0.6418117E+00 


0.6956997E+00 


4 


82 


0.3428497E+00 


0.3697114E+00 


4 


2841 


-0.4728359E+00 


0.1499717E+01 


5 


1 


0.1496636E-03 


0.4925571E-03 


5 


2077 


0.2547383E+00 


0.9045169E+00 


6 


1 


-0.1994387E-03 


0.4756932E-03 


6 


17 


0.1862375E-02 


0.8850131E-02 


7 


1 


0.3741456E-03 


0.4675257E-03 


7 


12 


0.1632629E-02 


0.9739530E-02 


8 


1 


-0.7509839E-04 


0.4525436E-03 


8 


44 


0.1126855E-02 


0.1116815E-01 


9 


1 


-0.3260709E-03 


0.4070134E-03 


9 


10 


0.1124477E-02 


0.7585276E-02 


10 


3 


0.1876750E-03 


0.3783614E-03 


10 


50 


-0.1077576E-02 


0.1420463E-01 



Table 1: Typical lump structures on a 10^ lattice, originally prepared at /3 = 2.3, after 
1000 constrained cooling sweeps. Left: action-based definition with a threshold t = 20% . 
Right: charge-based definition with a threshold t = 10% . 



We have used the standard notation F^^, := F^^t" and tT{tH'') = , where t"- (= Aa/2 
for 5^7(3)) are the generators of the su{Nc) Lie algebra. The action and the topological 
charge of a dilute gas of rij instantons and nj anti-instantons are given by 




For convenience, we also introduce the "normalised" action A by 

^ = T^E/^'^ (^r(^)^r(^) + Btix)B^ix) ) , (53) 

i a 

with the density 

= {E^{x)Et{x) + B^{x)Bt{x) ) (54) 
which gives rise to the global inequality 

A > \Q\ (55) 

and a local variant of it: 

— s{x) < q{x) < s{x) . (56) 

Let us briefly write down the lattice analog of the above quantities. Since this paper 
presents a first exploration of a new cooling procedure, we decided to use the most unam- 
bitious (unimproved) lattice observables. This is not doomed to failure since our cooled 
configurations are quite smooth. For example, there is no need for an improved definition 
of the field strength which uses other Wilson loops rather than the plaquette. In order 
to associate the field strength tensor with each lattice site, we introduce the clover-leaf 
average of the plaquettes attached to this site as indicated in the left panel of Figure O 
Each of the quantities Pi,, a represents an untraced plaquette as defined in (IT7I) with the 
difference that tadpole improved links 



Uf,{x)/uo , 
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Figure 6: Sample configuration obtained after 1000 constrained cooling sweeps applied to 
a Monte Carlo configuration originally prepared on an 16^ lattice at /9 = 2.4. Left: action 
density. Right: modulus of the topological charge density. 

were used in place of the standard links. For smooth configurations, i.e. if a? F^y{x) <^ 1, 
the clover average C^y(x) can be expanded with respect to the field strength: 

C^,{x) = 1 + ia^F;;,{x)e + ... . (57) 

It is the latter equation which we use to define the colour components of the field strength 
tensor on the lattice: 

a2F;,(x) = -22tr|rC^,(x)} . (58) 

As a first check, we investigated whether our cooled configurations still sustain the Ua{^) 
anomaly. To this aim, we calculated the topological susceptibility 

X = ^ {Q') (59) 

where the average is over configurations obtained after constrained cooling from the Monte 
Carlo ensemble. The result for 400 independent lattice configurations of size 16^ originating 
from a simulation at /3 = 2.4 is shown in the right panel of Figure [5] as a function of the 
cooling iteration number up to 1000 sweeps. We observe that the maximal value is reached 
after C(50) cooling sweeps. The asymptotic value is not essentially less. Using the string 
tension a = [440 Me V]^ to set the reference scale, our estimate for the asymptotic value 
(derived from the largest number of cooling steps) is 

X^/^ = (168 ± 3) MeV (60) 

is somewhat larger than the value reported in one of the earliest investigations [5] for the 
SU{2) gauge theory. There, minimising the action by a variant of the Metropolis algorithm 
(formally taking /3 — )■ oo) until the (first) plateau value of the action is reached, together 
with the twisted-plaquette definition of the topological charge density have allowed to 
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Figure 7: Left: action versus topological charge in a 2d-histogram (the colour-code gives 
the counts). Right: Fraction of the total action and of the total topological charge carried 
by the leading, action-based clusters up to ordinal number n (see the main text for a proper 
definition of the cluster construction and enumeration). 

estimate the topological charge of lattice configurations. Given the lattice sizes of the 
early simulations, we deem our result to be in agreement with the previous results for a 
SU{2) gauge theory. 

Figure [6] visualises the action density (left panel) and the modulus of the topological charge 
density (right panel) of a generic 16^ lattice configuration after 1000 sweeps of constrained 
cooling. We stress two important features: (i) both densities are well localised in the 
same, in space-time well separated objects; (ii) at first sight, the clusters seem to be (anti-) 
selfdual. We will quantify the selfduality, more precisely the degree of saturation of fl56l) . 
below. 

The left panel of Figure [7] shows a scatterplot of the total action versus the total topological 
charge. An ensemble of 4000 independent configurations prepared at /3 = 2.3 on a 10^ 
lattice has contributed to this result. We observe that, while the topological charge Q 
is quantised to a good extent as expected, the normalised action flS^ is widely spread 
in the scatter plot. This result is strikingly different from what a dilute gas pattern of 
instantons and antiinstantons would show: in that case, the action would be quantised as 
well (see ([52])). 

In the further analysis of the space-time structure, we have to identify the "lumps", i.e., 
clusters in a given configuration. We will employ different definitions of the clusters and 
compare the cluster properties arising in either case. 

Action-based clusters: 



22 



-6 



2000 



1500 



1000 



500 



a. I 



i 





4ump 



Figure 8: Left: histogram with respect to the summed topological charge of charge- 
based clusters. Right: histogram of g„ provided by the individual charge-based clusters. 
Note that the central peak, which reaches as high as 16000 counts was capped at 2000 
counts to fit into the figure. 



(i) At first, we scan the lattice for the maximal value Smax of the action density. Let this 
maximal value occur at point xo- 

(ii) Then, we label all points Xk which form a cluster connected to xq and with an action 
density s{xk) >tx Smax- All "connected" sites are connected to each other (directly 
or indirectly) by links. The threshold parameter t can be chosen < t < 1. The 
set of all these points Xk defines the first cluster Ci of the given configuration. This 
cluster gets the ordinal number n = 1. In later cycles of this cluster algorithm (see 
below) further clusters C„ with n = 2, 3, ... are defined. We then sum over the action 
density, a„ = '^^^c ■^(a;), and the topological charge density, = '^x&c Q'(^)' 
points which have been identified previously as part of the cluster Cn- Apart from 
action and topological charge, each cluster is characterised by Vn = "^xec 
number of points (cluster volume). 

(iii) We go back to step (i) and repeat the search for the maximal action density, but this 
time points which have previously been identified to be part of a cluster are ignored. 

After the cycle (i) to (iii) was repeated until all lattice points have been assigned to a 
cluster, we sort all the clusters according to their action a„ in descending order □ and store 
for the 10 biggest clusters a„, g„ and f„ for further analysis. 

^Therefore, these clusters are called action-based. 
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To give the reader an impression, a typical result from a single lattice configuration for a 
10^ lattice prepared at /3 = 2.3 before it underwent constrained cooling, is listed in the left 
half of Table [H Since the clusters are ordered and indexed according to the action, say, by 
n = 1 ... 10, we can ask how much of the total action and how much of the total topological 
charge is carried by a small number of leading clusters. To answer these questions, we define 
the cumulative quantities 

n 

i=l 
n 

Qn = ^Qi, (61) 
i=l 
n 

Vn = "^Vi ■ 

i=l 

Note that all three quantities implicitly depend on the threshold t parameter via the cluster 
definition. The right panel of Figure [7] shows that, with the above definition of the action- 
based clusters, a few of them already account for a good deal of both, action and topological 
charge of the whole configuration. 

In order to get further insights into the texture, we also defined the clusters using the 
topological charge density rather than the action density. Those clusters are obtained as 
follows. 

Charge-based clusters: 

(i) At first, we search for the maximal value of where q{x) is the topological 
charge density defined in fl^T]) . Assume that this maximal value occurs at point Xq. 
We then define the "sign" of the cluster to be defined first as a = sign(g(xo)). 

(ii) Then, we search all points Xk which are connected to Xq and which satisfy a q{x) > 
t |g(xo)|, where t is the threshold parameter. The set of all these points defines the 
first charge-based cluster Ci (with ordinal number n = 1) of the given configuration. 
In later cycles of this cluster algorithm (see below) further clusters Cn with n = 2, 3, ... 
are defined. As in the case of action-based clusters, the action = ^xgc^ ^i^)^ 
topological charge g„ = X]xgc„ volume Vn = J2xec„ ^ assigned to 
the cluster. 

(iii) We go back to step (i) and repeat the scan for the maximal value of \q{x)\ ignoring 
those points which already are part of a cluster. 

After the cycle (i) to (iii) was repeated until all lattice points have been assigned to a 
cluster, this time we sort all the clusters according to the modulus of the topological 
charge g„ in descending order and store the data of the first 10 clusters. 
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Figure 9: Scatter plots involving the volume Vn, the action a„ and the topological charge 
g„ of all individual charge-based clusters. 



A sample result of the cluster analysis with a threshold t = 10% applied to a single 
configuration on a 10^ lattice, Monte-Carlo generated at /3 = 2.3 and then constrained- 
cooled, can be found in the right half of Tabled! We then sum up the topological charges of 
the 10 leading charge-based clusters. By definition, summing up the charge of all clusters 
would yield the total topological charge Q of the lattice. We find, however, that - in 
particular if the threshold t is chosen smaller than 20% - the sum of the leading 10 lumps 
already gives a very good approximation to the total topological charge, i.e., Qio ~ Q- 
We emphasise the observation that applying a lower threshold parameter t leads to a more 
rapid convergence of the sum for the cumulative action and cumulative topological charge 
for the given configuration. 

The left panel of Figure [8] (left panel) shows the abundance of the cumulative topological 
charges Qsum = Qw for charged-based clusters using the minimal threshold t = for the 
cluster analysis of 10^ lattices prepared at /3 = 2.3. For this result, 4000 independent 
configurations have been analysed for their cluster structure. The lattice configurations 
underwent 1000 sweeps of constrained cooling. For small values of Qio, the cumulative 
topological charge peaks at integer values while for larger values a gradual shift to smaller, 
non-integer values is observed. We interpret these findings as follows: for small values of 
QiQ the lump texture is quite smooth allowing for an accurate estimate of the (integer) 
topological charge using the naive clover type of definition. For larger values of Qio, 
the configurations are more rough and call for an improved definition of the topological 
charge. With this caveat, connected clusters of action with charges up to Qiq < 6 have 
been observed. 

In order to investigate whether the topological charge Qium = Qn of the individual charge- 
based clusters is quantised, we have also calculated their abundance. The result is shown 
in the right panel of Figure [HI The cluster decomposition was performed using a threshold 
parameter t = 0. We observe a sharp peak around g„ ~ 0. This peak originates from a 
huge number of isolated micro-clusters (fluctuations) which contribute very little to the 
topological charge. More striking is the fact that the distribution is smooth for larger 
values of g„ with a maximum around ~ 2. No sign of an integer quantisation is visible 
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Figure 10: Left: the self duality measure R. Right: the inverse participation ratios and 
Z)"^. All quantities are shown as functions of the cooling sweep. 



as one would expect from a dilute gas of instantons. The charge-based clusters carry a 
broad spectrum of topological charges. Only the sum of their individual contributions is 
quantised. 

In Figure |9l we finally present scatter plots of data describing individual charge-based 
clusters with minimal threshold t = 0, involving their volume f„, action and topological 
charge g^. Each point in these figures corresponds to one cluster. There is a close correlation 
and monotonic relation between action and topological charge: the higher the action, the 
higher is the modulus of charge. This relation is not far from saturation, |g„| < |a„. 
Nevertheless, the space-time distribution of action and topological charge is not identical. 
Note that, in contrast to an instanton or meron gas configuration, the topological charge 
varies smoothly with the volume or the action of the cluster. We also observe a smooth 
correlation between the volume of a cluster and its charge and action, respectively. There 
are bounds an < lOvn/N'^ and < Ivn/N^- 

3.2 Selfduality and inverse participation ratios 

In the following, we will analyse to which extent the clusters are selfdual or anti-selfdual 
and which portion of space-time they do occupy. For this purpose, we introduce for each 
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144 


16^ 


18^ 


20^ 


R 


0.9244(1) 
2.1(1) X 10-3 
0.50(2) X 10-2 


0.9147(5) 
1.64(5) X 10-3 
0.36(1) X 10-2 


0.9379(5) 
1.18(3) X 10-3 
0.248(7) X 10-2 


0.9432(4) 
0.80(1) X 10-3 
0.160(3) X 10-2 



Table 2: The selfduality measure R and the inverse participation ratios for the field strength 
and the modulus of the topological charge density for 1000 constrained cooling sweeps. 



lattice configurations two normalised densities, "unit" vectors derived from s{x) and q{x): 

J^4U] = ^===, $^J-.J-. = l, (62) 
-^A^] = ^==' $^I).P. = l. (63) 
At first we notice that the ensemble expectation value 

'Y.J'MI^M) (64) 



vanishes upon the integration over the link fields. In order to find out whether a high 
field strength is linked to a high modulus of the topological charge density, we consider the 
expectation value R, 

^ < R = (Y. -^-t^] |^-[^]| ) - ^ ■ ^^^^ 

Note that R = 1 corresponds to a perfect uniform alignment between the modulus of the 
topological charge density and the field strength. Let us also point out that the overall scale 
of the colour electric and colour magnetic action density drops out from the definitions for 
R, J-' and V implying that not only regions of large action and charge density contribute. 
Also "small" , aligned fiuctuations can make a sizable contribution to these quantities. 

We have studied R for a 16^ lattice and for /3 = 2.4. The result is shown in the left 
panel of Figure [TU] as a function of the cooling sweep. We find that R gradually increases 
with cooling reaching values as large as 0.9344(6) after 1000 cooling sweeps. This finding 
indicates that the constraints, imposed upon cooling, do not hamper to reach a high degree 
of selfduality of the final cooled configuration. 

Furthermore, we studied the inverse participation ratios (IPRs) of the action density and 
the modulus of the topological density: 

'Y.^^) , ■■= (Y'^I) , (66) 
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Figure 11: Tlie inverse participation ratios {J^^) and {T>^) after 1000 constrained cooling 
sweeps as a function of the linear lattice size N {N^ being the lattice volume) for 400 
configurations of each size prepared a.t (3 = 2.4. 

If the configurations would be extremely localised, i.e., J-'xg would be not vanishing only 
at one site, we would find F*^ = 1. On the other hand, if J-' and V are smeared out over a 
sub-manifold of dimension d, we would obtain (with V = the lattice volume) 

Our result for F'^ and for a fixed number N is shown in the right panel of Figure [TOl 
We find that these quantities are increasing during cooling. This indicates an increase of 
the localisation of the configurations during cooling as expected. We also point out that 
the IPR for topological charge density, i.e., D^, is roughly twice as big as F^. This finding 
signals that the topological charge density is more localised than the action density. This 
could mean that the clusters of action possess one or more selfdual centres (giving rise to 
the localised topological charge density) and an halo of action made from more generic 
(not necessarily selfdual) fields. 

In order to study the effective dimension d of the cluster structure, the IPRs F'^ and -D^ 
(and also R) were calculated for configurations with the lattice sizes 14^, 16^, 18^ and 
20"^, all prepared at /3 = 2.4. The results obtained from an average over 400 independent 
configuration are summarised in Table [2] and in Figure [Til The results are fitted by 



1.9(1) X Ar-2-6(i)^ 



^ 13(1) X Ar-3-o(i). 



(67) 
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Note that not only the absolute values for the IPRs F"^ and D'^ are quite different as already 
noticed, but also the exponents are slightly different. At the same time, we observed a 
strong correlation between regions of large action density and of large topological charge 
density (brought to light by R). This indicates that there is a subtle difference in the 
arrangements of action density and topological charge density which are hardly visible by 
the naked eye (see Figure [6]). Let us also stress that the data fitted by fl67j) are obtained 
for a fixed number of cooling steps (i.e. 1000 sweeps). These figures might still change if 
the number of cooling sweeps is further increased. A systematic study of this effect is left 
to future work. 



4 Spontaneous chiral symmetry breaking 

Besides confinement, the spontaneous breakdown of chiral symmetry is another outstanding 
feature of QCD which shapes the hadron spectrum and therefore the structure of matter. 
The crucial test for spontaneous chiral symmetry breaking is a non-vanishing quark con- 
densate (qq). It is related to the eigenvalue spectrum p(A) of the fermion operator M[U] 
by virtue of the Banks-Casher relation: 

(qq) = 7rp(0) . (68) 

For simplicity, we adopt staggered fermions which avoid fermion doubling while keeping a 
subset of chiral symmetry intact. The massless Dirac operator is 

4 

M[U] = J2v,ix) [U,{x)6,+,,y - Ul{x-fi)5,.,,y] , (69) 

where the phase factors are given by 

r]^{x) = (_l)-i+-+-M-i . 

We have numerically estimated the spectrum of the low lying modes for configurations 
created on a 10^ lattice with /3 = 2.4. For this aim, we have calculated the 50 modes with 
A closest to zero, taking them as entries to a histogram. We then averaged the histograms 
over 400 independent configurations. The result for un-cooled configurations is shown in 
the upper left picture in Figure [121 As expected, the eigenvalue spectrum is non- vanishing 
at zero eigenvalues. We subsequently studied the eigenvalues spectrum after 10, 100 and 
1000 sweeps of constrained cooling. It turns out that cooling pushes the eigenvalues away 
from zero while a group of eigenvalues remains forming a peak near zero eigenvalue. It 
is striking that number of near-zero modes firstly decreases, but that the full strength of 
the spectral density near zero seems to be regained asymptotically. We finally compare 
the spectrum after 1000 sweeps of our new constrained cooling technique with that from 
standard cooling (see figure [T^ . While the near zero- modes are protected throughout the 
constrained cooling sweeps, the near-zero modes have practically disappeared after 1000 
sweeps of standard cooling. 
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Figure 12: The spectral density of the staggered Dirac operator as a function of its eigen- 
values for several numbers of constrained cooling sweeps. 
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13: The spectral density of the staggered Dirac operator as a function of its eigen- 
after 1000 sweeps of cooling. Left: for constrained cooling; right: for standard 



5 Conclusions 

Standard cooling finally approaches classical Yang-Mills configurations which lack impor- 
tant low energy features of the quantized theory such as confinement. As a tool for analysis, 
standard cooling has to be stopped by an ad hoc criterion. Although it might be possible 
by this procedure to generate an ensemble of configurations which grasp the essence of low 
energy Yang- Mills theory, the persuasiveness of the approach is fiawed since the properties 
of the configurations might still depend on the stopping parameters. 

Constrained cooling in the spirit firstly proposed in [13] and detailed in the present paper 
removes this deficiency. Here, as well as in [43], the respective physical constraint makes 
it unnecessary to monitor and stop the cooling process at all. Instead, the physical con- 
straints ensure that the configurations hardly change after many cooling sweeps and that, 
at the same time, the emerging (semi-) classical configurations retain the important IR 
properties. A natural choice for such a constraint is to retain the asymptotic form of the 
static quark- ant iquark potential. In this paper, we achieved this by demanding that every 
Polyakov line on the lattice, i.e., through each lattice site and in every direction, does not 
change during the cooling process. Despite of the huge amount of constraints, we find that 
smooth (semi-) classical configurations with very low action emerge (see the discussion in 
subsection I2.2l which quantifies this statement). Once the confining property of the (semi-) 
classical configurations is established, we find that other important low energy properties 
are conserved as well: the topological susceptibility is largely unaffected by constrained 
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cooling, and the investigation of the eigenvalue spectrum of the staggered Dirac operator 
indicates that these configurations provide spontaneous chiral symmetry breaking as well. 

A thorough analysis of the space-time structure of the configuration emerging from con- 
strained cooling has been performed. It turns out that these configurations can be described 
by connected clusters which are (anti-) selfdual to a good extent (see subsection 13.21 for 
quantitative details). We find that these objects are genuinely different from instantons: 
while the total topological charge of each configuration is quantised in integers (as it should 
be), the topological charge of individual clusters takes values which are continuously dis- 
tributed. Data showing the topological charge of individual clusters versus their action, 
see figure [9l seem to fall on a well-defined curve in the action-topological charge histogram 
(rather than populating isolated islands as it would be the case for instanton ensembles). 

To explore these many observables reported here with a high-level statistics, these first 
results have been obtained by using an extremely simple lattice set-up. In view of our find- 
ings, it is certainly worthwhile to further study the configurations exposed by constrained 
cooling employing an improved lattice action, an improved definition of the topological 
charge density and a more elaborate Dirac operator with maximal chiral symmetry. We 
also expect that extending the investigations to finite temperatures will provide valuable 
insight in the nature of the deconfinement phase transition and possibly in some features 
of the quark gluon plasma close to Tc. This, as well as an extension to SU{3), is planned 
for future work. 
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A Solving the constraint 

Using (1271) . the task is to determine the Lagrange multiplier such that the constraint in 
(B^ is satisfied. We firstly note that 

Af' = det{B,) + XMBlV,) + X% 
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and squaring the constraint equation (!25|) leads to a quadratic equation for A^: 

A fr(B^V ) + XK^.)-det(i?,)itr^(^,P,) _ 
A, + A,tr(S,P,) + 1 _ ltr\U,V,) " 

Since both U^, G SU(2), we point out out that 



1 



1 - -tT\U,V,) > 



We find two solutions: 



One of the solution is spurious, since we have squared the constraint equation to derive 
them. Indeed, calculating 



A/", . 1 



4 

rules out the solution A~ since J\f > 0. We also find the normalisation factor to be 
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